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Abstract. 

^^ . In this paper we combine the WMAP7 with lookback time and Chandra gas fraction 

Z) \ data to constrain the main cosmological parameters and the equation of state for the 

(■h ' dark energy. We find that the lookback time is a good measurement that can improve 

O ,. the determination of the equation of state for the dark energy with regard to other 



o 



o 



external data sets. We conclude that larger lookback time data set will further improve 
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5— i ■ our determination of the cosmological parameters. 
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£--~ ■ The seventh-year Wilkinson Microwave Anisotropy Probe (WMAP) data rigorously test 

the standard cosmological model placing constraints on its basic parameters. The 
WMAP measurements alone are not enough to break the degeneracy among some 
cosmological parameters or to place constraints on non-standard cosmological models. 
For example, measurements of cosmic microwave background (CMB) power spectrum 
alone do not strongly constrain the curvature of the universe characterized by the energy 
density parameter Qk- One needs to complement the CMB data with the luminosity 
or angular diameter distances measurements in order to constrain Q^ because the 
astrophysical distances depend also on the expansion history of the universe [1] . 

There are conclusive evidences that the universe is in a state of accelerated 
expansion. The Hubble diagram of Type la supernovae (SNela, [2j [3j H]), combined 
with CMB anisotropy measurements ([HE]), baryon acoustic oscillations (BAO) from 
the galaxy distribution data (0, [3, [9]), and galaxy cluster gas mass fraction 
measurements (pH, EH EG]) support the idea that we live in a spatially- flat universe 
where nonrelativistic matter make almost 30% of the critical energy density while the 
rest is an unknown component called dark energy, with negative effective pressure being 
responsible for the present phase of accelerated expansion of the universe |13j . 
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Constrains on dark energy density parameter VL^ and on its equation of state w (the 
ratio of pressure to energy density), can explain the nature of the repulsive force causing 
the acceleration of the universe. One possible explanation for this unknown component 
is an energy density constant in time and uniform in space. Such a cosmological constant 
(A) was originally postulated by Einstein to explain a static universe, later rejected when 
the expansion of the Universe was first detected and presently reinstated to account 
for the dark energy. Still, the computed value of A is expected to be 10 120 larger 
than the observed one. Another cosmological scenario considers that dark energy is a 
dynamical scalar field with a time varying equation of state. An alternative explanation 
of the accelerating expansion of the Universe is that general relativity or the standard 
cosmological model is incorrect. 

However, ground and space observations can not discriminate among different dark 
energy scenarios as the correct explanation of the observed accelerating universe: a 
cosmological constant, a dynamical scalar field or a modification of general relativity 

PUE]. 

The WMAP7-year data combined with other astrophysical measurements [lj place 
constraints on the dark energy. Assuming a flat universe (0& ~ 0), an accurate 
determination of the Hubble expansion rate {H ) helps in improving the limit of the 
equation of state of the dark energy [161 02] • m the paper of Komatsu et. all. [1] from 
the joint analysis of WMAP7+BAO+i?o i n the case of a time independent equation of 
state, a value of w = —1.10 ±0.14 at 68% CL was obtained. Furthermore, adding high-z 
supernova data to their analysis a more stringent limit was obtained, w = —0.98 ±0.053 
at 68% CL. However, this later result does not take into account the systematic errors in 
supernovae, which are comparable with the statistical errors [TBI fl~9] . Also, combining 
the cluster abundance and 5-year WMAP data, Vikhlinin et al. [20J found that for 
a flat universe w = —1.08 ± 0.18 at 68% CL. Furthermore, adding BAO [21] and the 
supernova data [22j, they found w = -0.991 ± 0.09 at 68% CL. 

In this paper, we perform a joint analysis of the CMB-WMAP7 data, constraints 
on Hubble expansion rate inferred from the age of astrophysical objects using lookback 
time method (LBT) and measurements of the gas mass fraction of relaxed clusters from 
Chandra X-ray observatory (Chandra). Our aim is to investigate for several parameters 
of a given cosmological model which combination of this three data sets puts better 
constrains. 

We choose the lookback time method because it has the advantage of using the ages 
of distant objects which are independent of each other so we can avoid biases present in 
techniques that use distances of primary or secondary indicators in the cosmic distance 
ladder method [T3] . 

Moreover, we use Chandra measurements because they currently provides one of 
the best constraints on Q m and have the advantage of being remarkably simple and 
robust in terms of its underlying assumptions [10] . 

The paper is structured as follows: in section 2 we briefly present the lookback 
time method, in section 3 Chandra gas fraction experiment, in section 4 we describe the 
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statistical analysis and present the data sets and finally, in the last sections, the results 
and conclusions. 

2. Lookback time method 



This time-based method uses ages of astrophysical objects (passively evolving galaxies 
or clusters of galaxies) which are independent of each other, to constrain cosmological 
parameters. The advantage of this method is that it avoids biases that are present in 
cosmic distance ladder method, offering an independent way to cross-check cosmological 
constraints obtained by using other methods [TJ. 

The lookback time is defined as the difference between the present age of the 
Universe (to) and its age at redshift z, t(z), 

r dz 

where H(z) is the Hubble expansion rate at redshift z. 

We use the spatially-flat cosmological standard model XCDM, X denoting the fact 
that the dark energy component has an time independent unknown equation of state. 
In this model, the Hubble expansion rate as a function of redshift can be written as: 



H(z) 



H [On(l + Zf + (1 - On)(l + zyV+°>Y'\ (2) 

where fl m is the matter energy density parameter and w is the equation of state for the 
dark energy. 



redshift formation 

GO Z f 

I 



redshift of the object today 
Zj z=0 

I I 



age of the universe 




Figure 1. The time evolution of the lookback time and the age of the universe 
as a function of redshift [23] (see also the text). 
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The observed lookback time t°^ s (zi) (see the diagram presented in Fig. 1), to an 
object i at redshift Zi is defined as: 

tt{zi)=tf°-t i {z i )-df. (3) 

where tQ bs is the measured age of the universe and ti(zi) is the age of the object (a 
passively evolving galaxy or a cluster) defined as the difference between the age of the 
Universe at redshift Zi and the age of the universe at the redshift of object formation, 

z f 

f Zf dz' 

U{ Zi ) = t{ Zi ) - t(z f ) = t L (z f ) - t L { Zi ) = j , 1 + z n H , z ,y ( 4 ) 

The third term in equation (3), df is a delay factor that encodes our ignorance 
regarding the formation redshift of the object Zf and is defined as: df = tQ hs — t^zp)- 

The lookback time method has the important feature that the age of distant objects 
are independent of each other and we can avoid biases present in techniques that use 
distances of primary or secondary indicators in the cosmic distance ladder method [7]. 

3. Chandra X-ray cluster gas fraction 

It is shown that the gas fraction(fgas) in X-ray luminous, dynamically relaxed clusters 
measured with the Chandra X-ray Observatory , helps in better constraining the 
cosmological parameters [101 128] . Following J25J EE], the matter content of the largest 
clusters of galaxies provides an almost fair sample of the total matter of the universe. 
Measurements of the ratio of baryonic and total mass in clusters of galaxies should 
closely match the ratio of Q^ and Q m . The matter energy density parameter fi m can 
be determined by combining measurements of the baryonic mass fraction in the largest 
galaxy clusters with accurate determinations of Q^h 2 from cosmic nucleosynthesis and 
CMB experiments, and a reliable measurement of the Hubble constant, H . 

The X-ray clusters gas fraction model fitted to a reference ACDM cosmology is [TO] 

KAjb(z) fn h 



- jACDM/VM L5 



l + s(z) \tt 



d A (z) 



(5) 



where d A {z) and <i^ CDM (2:) are the angular diameter distances to the clusters in the 
current test model XCDM and the reference standard cosmological model ACDM. 

dA = 7= sinh I \/Vtu / — — — I (6) 

H (l + z)^h V Jo H(z)J y ' 

The systematic uncertainties in the Chandra experiment are parameterized by the 

angular correction A, the non-thermal pressure support in the clusters 7, the baryonic 

mass fraction in stars s(z) = s (l + a s z), the bias factor b(z) = bo(l + ai,z) and an overall 

calibration parameter K for the residual uncertainty in the accuracy of the instrument 

calibration and X-ray modeling. 
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4. Analysis 

We perform our analysis in the framework of the extended XCDM cosmological model 
described by 7 + 1 free parameters: 

= ( Q b h 2 , Q DM h 2 , 9 S , t, w, n s , A s; df) (7) 

standard 

assuming uniform priors for all the parameters. 

Here Qbh 2 and QoMh 2 are the baryonic and dark matter energy density parameters, 
9 S is the ratio of the sound horizon distance to the angular diameter distance, r is the 
reionization optical depth, w = - is the equation of state parameter (p and p are the 
pressure and energy density of the dark energy), n s is the scalar spectral index of the 
primordial density perturbation power spectrum, A s is its amplitude at the pivot scale 
and df is the lookback time delay factor. Table 1 presents the parameters of our model, 
their fiducial values and the prior ranges adopted in the analysis. 



Table 1. The parameters of our model, their fiducial values and the prior ranges 
adopted in the analysis 



Parameter 


Fiducial value 


Prior range 


Q b h 2 


0.0223 


0.005^0.1 


VloMh 


0.105 


0.01^0.99 


S 


1.04 


0.5^ 10 


T 


0.09 


0.01 ^0.8 


W 


-1 


-2^0 


n s 


0.95 


0.5 ->■ 1.5 


ln[10 10 A s ] 


3 


2.7^4 


df(Gyr) 


1.5 


0^3 



We modified the CosmoMC Monte Carlo Markov Chain (MCMC) public package 
for our extended 9 parameter space to sample from the posterior distribution giving 
the following experimental datasets: 

The WMAP7 temperature and polarization CMB latest measurements pQ 
The Lookback time (LBT) measurements. From Ref. [21] we use the ages of 
six galaxy clusters in the redshift range 0.10 < z < 1.27. The standard deviation 
uncertainty for this age measurements is about 1 Gyr. We also use the age of 32 
passively evolving galaxies from Ref. [27] in the redshift interval 0.117 < z < 1.845. 
The error for this sample (la) is 12% of the age measurements. Therefore, we have 38 
measurements of U{zi) with uncorrelated uncertainties o~i. For the age of the Universe 
we use the WMAP estimate, t° bs = (13.69 ± 0.13) Gyr [5]. 
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Then, we compute the x 2 function for each cosmological model described by the 
set of parameters given in equation (7) as: 

2 m , ^ (t L (e)-tt% Zl ,df)) 2 (t (6)-cr , R x 

XlBtKV) = 2^ n 2 i -2 + ~2 > I 8 ) 

t=i * + <5 bs <5 bs 

where ti and to are the theoretical predicted values of the lookback time and of the 
age of the universe and t°^ s and t^ are the corresponding measured values. Also, o~i is 
refereing to the one standard error of the experimental data and cr° o bs is the uncertainty 
in the estimate of to- 

Chandra X-ray gas fraction. We used fgas data following the work from Ref. 
[10], [28] and [29]. We modified the CosmoMC Monte Carlo Markov Chain public 
package [30] to allow the use of fgas data when exploring the cosmological parameter 
space, implementing fgas module presented in Ref. [3Tj . 

Our implementation of this new module into CosmoMC package considers all 
seven parameters related to systematic uncertainties as being constant and having the 
following values: A = 0.2, 7 = 1.05, s = 0.16, a s = 0, b(z) = 0.824, a fe = and 
K — 1. This is different from the approach described in [10J, where either Gaussian or 
linear uncertainties are taken into consideration and after that, the seven parameters 
are added to the original list of 13 possible free parameters. In this manner, we manage 
to decrease the computational costs, especially when trying to constrain non-standard 
models using also the Chandra data. 

The x 2 function we used for fgas data has the form: 

2 , , _ f (f£ DM (z) ~ ft7s{z)f (r, - 0.214)' 
WJ "tT °7 9 l(z) 2 (0-022)2 {9} 

where f^ DM is the gas fraction fitted to a reference ACDM cosmological model and 
is given by (5) and fg%£(z), <J e /ga S (. z ) are t ne experimental values of gas fraction and the 
associated errors. In the above equation, r\ is the slope of the f gas { z ) i n a region with 
a radius for which the mean enclosed mass density is 2500 times the critical density 
of the universe at the redshift of the cluster, as measured for the reference ACDM 
standard cosmological model [10] . The difference between the experimental slope rj and 
its reference value is normalized to its expected standard deviation squared. 
WMAP7, LBT and Chandra data sets (WMAP7+LBT+Chandra) are combined by 
multiplying the likelihoods. We have performed a likelihood analysis using three 
cosmological data sets: CMB-WMAP7, lookback time and the X-ray cluster gas fraction. 

5. Results 

In order to see, for several parameters of the chosen cosmological model XCDM, which 
data provides better constrains we run the modified CosmoMC package on a parallel 
computing system by using 64 independent chains for the following combinations of data 
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sets: WMAP7, WMAP7+LBT, WMAP7+Chandra and WMAP7+LBT+Chandra. We 
impose for each case the Gelman & Rubin convergence criterion [32]. 

We present in Table 2, the mean value and la error at a 68% CL, obtained from 
our analysis, for the following cosmological parameters: Vt^h 2 , Qrjuh 2 , fi\, fi m , n s , w, 
Ho, age(Gyr), log[10 10 A s ]. In the first column, we present the results for the reference 
WMAP-WCDM modeljj and in the second column, our simulation with WMAP7 data 
alone. Then, in the next two columns we present the results obtain from the combination 
of WMAP7 with LBT and Chandra data and finally, in the last column, we show 
the results when all three data sets are joined. Comparing with the WMAP-WCDM 
reference model, we conclude that our computation best constrain the finish 2 , fi\, Vt m 
parameters as expected, if we use all three data sets (WMAP7, LBT, Chandra) and the 
parameters Qbh 2 , n s , w, Hq, Age(Gyr), log[10 10 A s ] are better constrained with WMAP7 
and Chandra data. 

In Fig. 2 we show the correlation between some cosmological parameters, 
(Q m ,Q A ), (Q m ,w), (H ,Q m ), (H ,w), (Age/Gyr,w), (Q\,w). This are the joint 
two-dimensional marginalized distributions with 68%, 95%, 99% CL for the following 
combinations of data sets: WMAP7 data only, WMAP7+LBT, WMAP7+Chandra and 
WMAP7+LBT+Chandra. The plane Q m — Q A is better constrain when we combine all 
three data sets. Also for others correlations, the WMAP7 and Chandra combination 
is the optimal one. Although adding LBT data set to the analysis we obtain better 
constraints upon all cosmological parameters, the central values of the distributions 
are shifted when comparing with WMAP-WCDM reference model, so we prefer to use 
only the WMAP 7+ Chandra data sets when we want to see the correlation between the 
equation of state of dark energy and other cosmological parameters: Q m , H , age/Gyr, 
Q\ and also, the correlation between H and Vt m . 

We have shown that, for some parameters of a cosmological model with an unknown 
time independent equation of state, adding LBT to Chandra and WMAP7 data, leads 
to an improvement of their error bars. In the future, we expect that a larger LBT data 
set will further improve the constrains on cosmological parameters. 



| http://lambda.gsfc. nasa.gov/product /map/dr4/params/wcdm^z_lens_wmap7.cfm 



Table 2. The mean and lc uncertainty for cosmological parameters in WMAP7, LBT and Chandra combined data analysis 
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WCDM-WMAP7 


W7 


W7+LBT 


W7+CHANDRA 


W7+LBT+CHANDRA 




Q b h 2 


n (122 P.S+000063 


n 099 ^q+o- 00058 

U.UZZOy_ .00058 


n 091 07+0. 00073 
U.UZIO/ -0.00068 


02202+° 00047 


0.01909+°;™ 
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Figure 2. Joint 2D marginalized distributions at 68%, 95%, 99% CL showing 
the correlation between (fl m ,fl\), (Q, m ,w), (Ho,Q m ), (Hq,w), (age,w), 
(CIa,w). The blue contour is for WMAP only data set, the red contour for 
WMAP+LBT, the yellow for WMAP7+Chandra and the green contour for 
WMAP7+LBT+Chandra data sets. 



6. Conclusions 



The aim of our paper was to analyze how the LBT and Chandra data sets combined 
with WMAP7 measurements can improve the determination of the equation of state for 
the dark energy and also of the other parameters of the standard cosmological model. 
We choose the LBT external data set because it contains information independent of 
each other and we can avoid biases present in the cosmic distance ladder method and the 
Chandra data set because it provides one of the best constrains on the matter density 
parameter Vt m . 

Our contribution consist in implementing LBT and Chandra new modules into the 
public Monte Carlo Markov Chain package. We run the modified CosmoMC package for 
the following combinations of data sets: WMAP7, WMAP7+LBT, WMAP7+Chandra, 
WMAP7+LBT+Chandra. We found that, the physical baryonic density parameter, the 
scalar spectral index, the equation of state for the dark energy, the Hubble expansion 
rate, the age of the universe and the amplitude at the pivot scale are best constrain 
using WMAP 7+ Chandra data sets. All the other considered parameters, the physical 
dark matter, dark energy and matter density parameters and the plane Q m — fix are 
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best constrain when we combine WMAP7+LBT+Chandra. 

Moreover, we conclude that the looking back time is a trustful measurement and in 
the future a larger LBT data set can further improve the determination of parameters 
we investigate. 
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